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Suiuinary 


Wc  present  a  new  formulation  of  the  Taylor  series  for  a  class 
of  functions  which  can  be  expressed  as  the  product  of  a  purely 
angular  and  purely  radial  part,  viz.. 


GCrJ  =  Y^^Cr)f(r) 


in  which  is  the  spherical  harmonic  function.  The  formulation 

is  based  upon  the  use  of  the  Fourier  transform  of  the  function  to 
be  expanded.  The  general  expression  for  the  Taylor  series  which  we 
obtain  is  similar  in  form  to  the  Laplace  expansion  for  the  Coulomb 
potential.  Thus,  Taylor  scries  can  be  developed  for  arbitrary 
functions  along  lines  which  arc  similar  to  the  multipolar  expansions 
in  tlie  electromagnetic  theory.  Such  expansions  are  easily  adapted 
to  the  symmetry  of  a  collection  of  sources. 


1.  Introduction 


dur  purpose  in  this  paper  is  to  present  a  new  and  nowerrul 
method  of  handling  the  expansion  of  a  class  of  functions  in  a  Tavlor 
series.  The  elements  of  this  class  arc  simply  those  functions  which 
can  iic  expressed  as  a  product  of  a  purely  radial  and  purely  angular 
part :  viz.  , 


(^■'1  >■' (1-1 

in  whicli  i^pherical  harmonic  function.  We  jiroceeil 

to  derive  the  general  term  of  the  I'aylor  series  to  any  arbitrary 
order.  We  arc  able  to  show  that  given  the  Fourier  transform  of 
any  wel  1 -betiaved  function  (or  generalized  function"),  it  is  possible 
to  obtain  the  specific  form  of  an  arbitrary  term  in  the  series. 

Tliese  results  extend  and  generalize  an  initial  effort  whicli 
was  presented  elsewhere  by  us  (Schmidt,  Pons,  and  McKinley,  198(1"). 

In  particular,  by  more  laliorious  means,  we  derived  specific  formulae 
for  the  first  and  second  order  terms  of  the  Taylor  scries.  The 
method  which  we  used  in  that  earlier  work  we  found  not  to  generalize 
easily  to  the  expression  of  a  general  term  for  the  Tai-lor  scries. 

The  results  of  this  paper  contain  exactly  the  results  reported 
previously  (Schmidt,  ct  al . ,  1980). 

We  find  several  reasons  for  undertaking  the  development  of  this 
new  expression  of  the  Taylor  series.  Foremost  is  our  need  to 
consider  the  displacement  of  a  particle  about  a  point  of  ecui  i  1  i  I' r  i  urn 
which  simultaneously  is  also  a  centre  of  s>’mmetr\'  for  a  collection 
of  .sources  (Schmidt,  et  al.,  1980).  Thus,  we  need  to  consider 
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a  form  of  the  Taylor  series  expansion  of  an  arbitrary  function 
which  exploits  the  symmetries  of  the  system  in  the  same  manner 
as  is  customarily  clone  for  the  electrostatic  potential  through 
the  use  of  the  Laplace  expansion  (cf Jackson ,  19h2).  These 
symmetries  are  most  effectively  handled  through  the  use  of  expansions 
which  depend  upon  the  spherical  harmonic  functions.  Hence,  it  is 
natural  to  consider  the  use  of  Fourier  transforms. 

Often,  and  for  molecular  physics  in  particular,  it  is  necessary 
to  consider  approximations  of  complicated  functions.  This  is 
especially  true  of  potential  energy  functions  which  arc  derived 
from  molecular  integrals.  The  analysis  of  the  mechanics  of  molecular 
vibrations,  for  example,  commonly  depends  upon  approximations 
which  replace  the  exact  potential  energy  functions  with  harmonic 
replicas.  It  is  universally  known,  of  course,  that  these  harmonic 
approximations  are  adequate  only  when  applied  to  the  analyses 
of  the  molecular  low-lying  vibrational  states. 

It  is  frequently  of  interest  to  examine  higher  order  terms  in 
order  to  assess  the  accuracy  of  the  harmonic  approximation.  In 
some  cases,  higher  order  terms  are  required  in  order  to  establish 
limits  of  stability  for  mechanical  systems.  For  a  relatively  simple, 
complete  function,  it  is  possible  to  extract  cartesian  harmonic  and 
fiigher  order  terms  as  individual  terms  in  a  Taylor  series.  Realistic 
potential  energy  functions,  however,  often  contain  fairly  complicated 
angular  dependencies.  These  angular  dependencies  reflect  tlie 
complicated  character  of  the  environment  which  surrounds  a  particular 
particle  (atom  or  molecule)  of  interest  [see,  for  example,  Briels 
(1980)1.  Pot"  these  functions,  the  forming  of  the  Taylor  series 
ceases  to  he  a  simple  matter. 
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The  usual  representation  of  the  Taylor  series  presents  the 
determination  of  the  coefficients  as  a  sequence  of  differentiations. 
These  operations  are  applied  to  the  function  the  expansion  of  which 
is  required.  It  is  possible  in  principle  to  make  repeated  use  of 
Rose's  (19S7)  formula  for  the  gradient  to  determine  these  coefficients 
in  terms  of  the  basis  vectors  of  the  spherical  tensors.  In  reality, 
however,  any  actual  attempt  to  apply  that  formula,  even  to  the 
simplest  of  functions,  quickly  becomes  unmanageable.  It  is  tedious 
even  to  determine  the  harmonic  terms  for  the  series. 

In  contrast,  we  are  able  to  obtain  the  coefficients  of  the  Taylor 
series  after  the  evaluation  of  only  a  single  form  of  radial  integral. 
The  strength  of  the  method,  therefore,  is  a  transparent  compactness 
and  flexibility  in  terms  of  ease  of  application.  This  strength  is 
illustrated  by  several  examples  in  a  separate  paper  (McKinley  and 
Schmidt ,  198_  b) . 

The  outline  of  this  paper  is  the  following.  In  section  2 
we  consider  the  single  centre  expansion  of  a  function  into  a  Taylor 
series.  A  particular  form  of  our  result  replaces  differentiations 
with  integrations.  For  complicated  functions,  which  possess  well- 
defined  Fourier  transforms,  these  integrals  prove  easier  to  consider 
and  evaluate  than  is  the  case  of  the  consideration  of  the  multivariate 
differentiation  of  the  function  directly. 

In  section  ,  based  on  the  use  of  the  results  of  section  2,  we 
present  a  derivation  of  Rose's  (1957)  formula  for  the  gradient 
operation.  A  consistency  between  Rose's  work  and  ours  is  established. 
In  addition,  we  use  the  differential  raising  and  lowering  operators 
for  the  spherical  Resscl  functions  to  re-establish  a  differential 
form  for  the  laylor  series  coefficients.  This  differential  form 
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is  consistent  with  the  angular  dependencies  of  tlie  integral  form 
and  with  the  svmmetries  of  the  spherical  harmonic  functions.  It 
is  clear  from  the  differential  form  obtained  that  in  many  instances 
the  integral  representation  is  indeed  the  simpler  way  to  approach 
the  determination  of  coefficients  in  a  Taylor  series. 

In  section  4  we  consider  the  Taylor  series  as  an  expansion 
about  two  centres.  It  is  possible  to  consider  the  general  develoji- 
ment  of  any  well-behaved  function  as  an  expansion  about  two  or 
more  centres.  The  Carlson-Rushbrooke  (1950)  expansion  of  the  free- 
space  Green  function  is  such  an  example.  In  section  4  we  establish 
a  parallel  treatment  of  the  Taylor  series.  We  limit  our  considerations, 
however,  only  to  two  centres.  In  this  section  we  also  show  that  as 
tlie  two  centres  coalesce,  the  resulting  form  is  consistent  witli  the 
single  centre  form  of  the  Taylor  series. 

Finally,  in  the  last  section  we  mention  the  limiting  forms 
for  scalar  functions.  Whereas  for  the  general  case,  for  a  general 
angular  dependence,  the  integral  form  of  the  Taylor  series  seems 
to  be  the  simpler  route  to  follow,  for  a  scalar  function,  we  find 
the  resulting  differential  form  is  equally  simple, if  not  simpler,  to 
use.  Our  form  of  the  Taylor  series  in  its  differential  representation 
for  a  scalar  function,  however,  does  not  bear  a  close  resemlilanco 
to  the  cartesian  differential  form.  Ou’"  form  preserves  the  angular 
deiiendenc  ies  as  arguments  of  the  Legendre  polynomials.  nifferentia 
tion  is  ajiplied  only  to  the  radial  part  of  the  function. 
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.  L'hc  Taylor  Scr i cs  for  the  Expansion  about  one  Centre 


We  develop  the  Taylor  series  for  functions  \vhich  arc  separable 
into  purely  radial  and  purely  angular  parts,  viz.. 


G(r)  =  F(r)Y^^(^) 


(2.1) 


where  Yj^^(r)  is  the  spherical  harmonic  function.  Still  more  general 
functions  could  be  expressed  as  superpositions  of  such  functions. 

The  function  G(r)  admits  a  Fourier  transform: 


Girl  =  (2  it) 


-  .5 


d'k 


f (k) Y^^ (k)cxp(- ik- 


r) 


where  f(k)  is  given  by 


f  (k) 


4  IT  i 


dr  r^F(r)j^(kr) 

0 


and  J^(x)  is  the  spherical  Bessel  function  of  the  first  kind. 

The  function  at  the  displaced  point  r  +  c  is  given  by  the 
vectorial  Taylor  scries  (Band,  1959) 


‘■ir  +  c)  =  I  i  fc-V)"G(r)  .  (2.41 

n  =  0  ''  ' 


Consider  now  the  general  term  in  this  series  with  G(r)  given  by 
its  t  .  : 
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■tc*V)  c:(,r) 


1 


(2Tr)  ! 


d^k 


r  (.k)  (k1  ( -  i  c  •  k)  ’’exp  [  -  i  k  •  r  ) 


j3|.  i^nfCklV,  (k)  (c-kl'Vxpl- ik-rl  . 

C2Ti)^n!J  -  - 


Here,  c*k  is  the  cosine  of  the  angle  hetween  the  unit  vectors  c 
The  nth  power  of  this  quantity  can  be  expanded  in  l.egendre  poly 
nomials  of  the  same  argument  {.Morse  and  I'eshback,  19S3) 


/N  /N  /\  /N 

(c-k)"  = 


The  coefficients  arc  given  by 


A 


2^+1 


nil 


r+1 


-1 


dx  X  l\(x) 


0,  l>n  or  n-H  odd 


(.2£^-l)n!  Cn-il-*-])  !  ! 
(n-il  +  i)  !  (n  +  Jl  +  l)  !  ! 


,  !l<n  and  n-il  even 


We  now  introduce  the  Rayleigh  expansion, 


cxp('ik*rl  =  { -  i )  ^' {  2b+ 1 )  j  .  (kr )  1>  (k  •  r )  , 


and  the  addition  theorem  for  the  spherical  liarmonic  functions, 


21, +  1  1.  ’ 

M 


2.  SI 

a  lul  k  . 


:.h) 


(2.71 


12.8) 


so  that  we  have 


-s- 


cxpi-ik'r)  -  U  y  (  -  1  )  i  I  (  kr  t  V|  I  V|  r 
L  ,  M 


?  ,m 


(2.11) 


With  these  substitutions,  Equation  (2.5)  becomes 


1  n  ■’  (' 

2r(>;-'?)  <Url  «  4- 


L  ,M  ,Jt  ,m 


*  ^  * 


X  dk  k-  dil^k“f{k)  j^(kr)Y^^(k)Yj^,(k)Y5^,Jkl 


(2.12) 


The  angular  integral  is  the  well-known  integral  over  three  spherical 
harmonic  functions  (Rose,  1957) 


i  ^•■-i.^„(k)Y;'„(k)Yt.,.(k)  = 


J  k  Au  I.M  £m 

In 


(l.niO  1  XO)  (l.f.Mml  >u) 


(2.15) 


where  (I,!Mm|Aii)  is  the  El ebsch - Gordan  coefficient.  We  also 

intro'’uce  a  quantity  defined  by  the  remaining  radial  integral: 


1.1  (r)  =  -  dk  k"^^f(k)j,  (kr) 

(27t)^Jo 


(2.14) 


Altogether,  we  have  found  the  following  expression  for  the  general 
term  in  the  Taylor  series: 


-(c-V)'‘n(r)  =  (47t) 


.5/2  c" 


1 


I.  ,Mji  ,m 
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2L+1 


(2£+1) (2X+11 


1/2 


(LILOO  I  XO)  (L£Mm|  Xy)  (r)  . 


(2.15) 


The  most  frequently  used  cases  will  he  those  of  the  lowest 
order.  For  n=l,  il  can  only  be  1  and  Aj^  =  l.  We  have  then 


c-VGCr)  =  j 

L  ,  M ,  m  ^ 


1/2 


(LlOO |X0) (LlMml Xu) Ij j  (r) . 


(2.16) 


A  special  application  of  this  is  to  find  an  equilibrium  condition, 
which  can  be  expressed  as  C‘VG(r)  =  (1,  for  any  c,  or  as 


0  =  J  (-i)W,^.(r)(2L+l)^/^(L100|X0)(LlMm|Xu)I,,  (r).  (2.17) 

L,M 


For  n=2,  H  can  be  either  0  or  2  with  A2q=1/3  and  = 


We  have 


then 


0  0 


T  I  (■i)‘'*b,„(r)Y,„(c) 


L  ,M,m 


■2m' 


2L  +  1 


5(2X+1)J 


xl2L(r) 


(c)l2;^(r) 

1  /2 

(L200 ( XO) (L2Mm| Xu) 

(2.18) 


These  lowest  order  terms  are  obviously  easy  to  evaluate,  hut  even 
the  general  term  [F.quation  (2.15)]  is  not  intrinsically  more 
difficult.  This  simplicity  traces  back  to  Equation  (2.11)  which  is 
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an  especially  simple  example  of  the  expansion  of  an  invariant  into 
irreducible  forms  (I'ano  and  Racah,  1959). 

3.  The  Formula  for  the  Ciradient 


In  this  section  we  verify  that  the  first  order  term,  eqn  fZ.lh), 
is  eciuivalent  to  the  usual  gradient  formula  (Rose,  1957). 

We  begin  by  evaluating  the  Cllebsch-Gordan  coefficient  (blOOl'O) 
in  expression  (2.1b).  This  coefficient  restricts  the  parity  of  I. 
as  well  as  limits  its  range.  We  have 


1/2 

(A-1,100|A0)  = 


(A  +  1,  J00|  AO) 


[iattJ  ' 


(.3.1) 


all  other  values  of  L  yield  zero.  We  now  write 


c*VG(r)  =  (-i)^^^^^ 


3/2 


I 


2AVT 


1/: 


/3  M,m 


(A-l,]Mm|Ay)Y^_j^j(r)Yj^Jc) 


(3.2) 


We  handle  the  angular  and  radial  factors  separately.  The  standard 
definitions  of  the  vector  spherical  harmonic  functions  are  (rose, 
1957) 


M  ,m 


(3.3) 


where  the  quantities  are  the  basis  vectors  for  the  spherical 
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tensors  (Rose,  1957).  The  scalar  product  of  (7.3)  with  the 
vector  c  i s 


J  UtUM«|X„)Y,,,„,(hc-„, 

M  ,m 


=  ./4¥77  I  (AillMmI  An)Y^^j^,(^)Y^^^^(c)  . 
M ,  m 


(3.4) 


On  the  other  hand,  for  the  radial  part  we  consider  the 
following.  In  the  definition  of  standard 

operators  for  lowering  and  raising  the  index  of  tlie  spherical 
Resscl  functions  (Morse  and  Feshhack,  1953,  and  Infeld  and  Hull, 
1951).  Thus, 


<lk  (kr) 

) 

o 

dk  k*f(k) 


^  0 


fA+1  ^  d 

[~r  HT 


*  TOFT  h"'’'’ 


dk  k^f(k).i^(kr). 


(.3.5) 


Similarly, 


(2TT)Mj^^^(r)  =  I  dk  k^f(k).i^^3  (kr) 


A  _  d 
r  cTr 


dk  k^  f (k) (kr) . 


(3.(0 


At  this  point,  we  identify  the  inverse  of  eqn  (2,5) 


F(r) 


_  (-i) 


271  ^ 


dk  k^f(k)j^(kr) 


(  3 .  "  ) 
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llpon  sul'st  i  tut  ion  of  all  of  these  expressions  into  C(|n  (3.2), 
we  find  exactly  the  standard  formula  for  the  j’radicnt  operation 
(Rose,  1057); 


This  new  form  in  terms  of  Touricr  transforms  could  have  been 
anticipated,  because  the  radial  operations  in  cqn  (3.8)  are  uniquel)- 
those  associated  with  the  splierical  Bessel  functions. 

bqn  (3.8)  completes  the  connection  with  Rose's  formula  for 
the  (iradient.  We  now  extend  the  analysis  piven  above  to  the 
consideration  of  an  arbitrary  term  in  the  Taylor  series.  Thus,  we 
consider  the  transformation  of  an  arbitrary  term  from  an  intepral 
to  a  differential  form  in  the  radial  component  only. 

For  a  piven  function  (Ifr),  A  is  s))ecified.  J'or  a  piven  value 
of  the  order  n,  the  values  of  I  also  are  specified.  Tlie  ('lebsch- 
r.ordan  coefficients  fix  the  parity  and  ranpe  of  the  f-values. 

Thus,  it  is  possible  to  see  that 

I,  =  A  +  q  (3.',') 


with  the  values  of  c)  fixed  and  limited  In-  the  conditions  mentioned. 
The  summation  over  L  in  the  expression  for  an  arbitrary  term  in  the 
Taylor  series  is  rewritten  as 
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(c-Vl”('.(r1  =  1471) 


3/2  c' 


I  (-1)"“— 


II ! 


n ! 


lie 


e  ,m 


[  (2e+n  (Ze  +  I  )  1 


,TT77  e.iu 


x  I  I 

q  =  0  M 


(-i)^V2(’X+q")+T(A+qeOO|  AO)  (A+qeMm  |  Aq)  ^  f  r) 

A 

+  (-i)'''/2(A-q)VT(A-qeOOl  AO)(A-qeMm|  Aq)Y^_^^^,(r)rj^^;^.^(r) 


(3.10) 


In  fact,  tlic  summations  run  over  q=0,2,4,6...  or  q  =  l,5,5,... 
depending  upon  tlie  values  of  n  and  A.  The  factor  togetlicr  \\ith 

the  Cl cbsch-flordan  coefficient  automatically  sort  out  the  terms 
which  survive  in  the  summation.  Hence,  there  is  no  need  to  display 
any  specific  form  of  the  sorting  process. 

We  now  use  the  following  lowering  and  raising  operators 


=  (-1) 


,A  +  q 


rdr 


^  -Ar,-A. 

r  {k  j 


(kr)} 


(3.11) 


and 


q- A-1 


d 


■^^^{k'^j^(kr)} 


(3.12) 


to  write 


d  -Af 


dk  k"''^^"f(k)j^(kr)  (3.13 


and 


-A-1  f  d  VLa  +  1 
rdr 


dk  k"''’*"f(k)j^(kr)  (3.14) 
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For  a  given  value  of  q  in  1  ,  (rl  we  have  adiusted  tlie 

n  ,  A  +  q 

index  on  the  Bessel  function  to  a  value  of  X.  It  remains  tlierer(Mc 
onlv  to  lower  tl)e  magnitude  of  the  exponent  of  k  in  the  integraiul 
from  to  k“.  This  we  do  as  follows.  First,  note  that  from 

the  various  conditions  on  the  parities  of  the  indices,  n-q  will 
always  be  even.  Thus,  we  can  reduce  n-q  to  a  zero  value  in  steps  of 
2.  From  the  eigenvalue  equation  which  the  Bessel  function  j 
satisfies,  we  write 


X^i  (x) 


(3.  IS) 


The  right  hand  side  of  this  equation  defines  an  operator  which  lowers 
tiie  exponent  on  k  (x  =  kr)  by  2.  The  use  of  this  operator  (n-q)/2 
times  yields  the  desired  result: 


dk  k"'''‘'“f  (k)  j ,  (kr)  = 


n-u 


r" 


dk  k'f(k)i^ (krl 


2tT  ^  i 


.  X 


n-q 

'•1  /-  „  2  ti  ^  1  2 


UX  +  1)  -  r'i'^inr). 


(3.1b) 


The  final  expression  for  the  arbitrary  term  in  the  Taylor 
series  is 


nt 


i  (c.V)"c(r)  =  /4it  ^  I  (-i)"  - 

~  ~  f,,m  (  (2X  +  1  (2J.  +  1)  1 


T7I 


'  I  I  i‘' 

q=0  M 


/2(X+q)+l  (X+q?,00  j  XO)  (X  +  qR,Mm|  ^  F )  Y ^ ( r )  r 


X+(j 


d 

rd  r 


-A 

r 


+  ¥ 


^rx-q)  +  l(i-qf,n()|XO) 


(3.1^) 


1.  Hic  r.xpans  ion  ol  a  I'unction  in  a  Tavlor  Series  about  two  Centres 


lor  some  apj)  1  i  cat  ions  it  is  useful  to  consider  the  function  of 
ei|n  tJ.lJ  to  depend  upon  two  points  ri  and  ra  ,  either  of  which 
may  he  displaced  independently.  We  set 


r  =  r,  -  !•; 


(4.1) 


and  redefine 


(i(ri  .ra  )  =  r)l'(r)  . 


(4..:) 


I'liis  function  still  possesses  a  i-ourier  transform  as  exjiressed  hy 
eqn  (2.2),  (2.3)  and  (3.7J.  Before  we  evaluate  this  function  at 
displaced  points,  however,  we  present  an  alternative  Fourier  transform 
which  emjdiasizes  the  two-centre  character. 

Into  the  expression  for  the  Fourier  transform 


(:(ri,r2)  =  (2ti)  d’k  f  (k)  Y.  ( k)cxp  ( -  i  k  •  r  i  )  exp  ( ik*  r2 ) 


(4.3) 


we  substitute  Luw  Rayleigh  expansions  of  the  form  (2.10).  The 
angular  integrations  over  the  three  spherical  harmonic  functions 
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G(.ri  ,r2) 


^f. 


dkk' 


dQ,fCk)Y  (M  t 

4ti  Li  ,Mi  ,  L2  ,M2 


3/7 

1T-"^“  LiMiL2M2 
1/2 


^  (LiL200lA0)(L,L2MiM2|Xy) 


.00 

cj  dk  k^  f  (k)  j  (kri  )  (kr2 ) .  (4,4J 


In  order  for  eqn  (4.4)  to  be  consistent  with  cqn  (.^.7)  and  (4.2) 
for  arbitrary  functions,  we  require  the  following  double  addition 
formula  for  spherical  harmonic  and  Bessel  functions: 


1/2  V  .X+L2-L, 


Y^^(r)j^(kr)  =  (4tt)^'"  I  1 


(L1L2OOIXO)  (L1L2M1M2  |Xvi) 


LiMiL2M2 


(2Li+l) (2L2+1)1 

zrn 


1/2 


A  similar  form  for  real  spherical  harmonic  functions  has  been 
derived  by  .Johnson  (1973).  It  is  possible  to  show  easily  that  the 
right  hand  side  of  (4.5)  reduces  to  the  form  of  the  left  hand  side 
when  r2=0  and  r=ri.  Additionally,  when  X=y=0,  the  right  hand  side 
of  (4.5)  can  be  shown  to  reduce  to  Yo 0 (r) j 0  ikr)  with  the  use  ol  the 
appropriate  Clcbsch-CIordan  coefficients  and  the  addition  theorem 
for  the  spherical  h.irmonic  functions.  Thus,  the  two-centre  expansion 
is  consistent  with  the  more  familiar  single-centre  forms. 

We  proceed  to  evaluate  the  two-centre  function  at  the  displaced 
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points  rj+a  and  rj+b.  We  have  as  the  extended  vectorial  Taylor 
series,  the  expansion 


GCri+a,r2+b)  =  I  i(a*Vi  +  b  •  V2 ) (r  i  ,  r  2 )  . 

'  '  '  ~  n=0  ~  ~  ^  ~  ~  ' 


( 4  .  b  ) 


The  general  term  in  this  series  is  to  be  evaluated  with  G(r],r2) 
expressed  by  its  Fourier  transform  [eqn  (4.3J].  We  have 


-^(a*Vi  +  b*V2)  G(ri,r2) 

■N.'v  -v. 


=  (2tt)'^| 


d'k  f(k)Y^^(k) 


X  p^[-ia*k  +  ib-k]*^  exp  ( - ik*r  1 )  exp  ( ik •  r  2 )  . 


(4.7) 


In  this  expression  wo  use  the  binomial  expansion  and  the  expansion 
in  irreducible  forms  (2.11)  to  replace  the  quantity  in  the  brackets: 


[a-k  -  b-kl"  =  (a-k)“''<(-b.k)'1 


0  I  ^q^^  ^  ^  (Z^i  +  l)  (2J12  +  1  )■ 

q  ,  ii,  1  , 2  m  1  ,  m  2  '  '• 


"'^'n-q.X./q  ,«,2'^«-inii  R  i  mi  P,2m2  ?  2mn  ’ 


(  1  .  .S  ) 


The  coupling  rule  for  the  spherical  harmonic  functions  is  now  used 
to  reduce  the  number  of  these  functions  in  the  product  (Rose,  1F571: 


fs  ^  ft  ^ 

4iTTmT~^J 


><(^1  e2mim2  I  . 


(  4  .  P  ) 
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The  substitution  of  two  Rayleigh  expansions,  and  the  use  again  of  the 
coupling  rule  together  with  eqn  (4.8)  and  (4.9)  leads  to 


1  [a-Vi  +  b-V2]''G(ri,r2)  =  I 

q,£i  ,%2  ,nii  ,in2  ,Li  ,L2  ,Mi  ,M2  ,L,M 

1/2 

''^L2M2  *'*^2^  ,  (2Jii+l)  (2jiJ  +  1)  (2£+'lfT2L+l))  ( «- 1 «- 2  00  |  £0)  (i.  i  J.  2m  i  m2  |  i'.m) 


x(LiL200|L0)  (L1L2M1M2  |LM)|  dk  (krj)Jj  ^  (kr2) 


4ti 


(4.10) 


The  angular  integral  is  evaluated  as  before  (cf.,  eqn  (2.13)),  and 
we  introduce  another  quantity  defined  by  the  remaining  radial  integral: 


-00 

J„,  ,  (ri.rz)  =  — ^ —  dk  k"'^^f(k)j,  (kri)j  (kr2). 

nLiLz  ~  ~  (-2^)3  Li  L2 


(4.11) 


Thus,  we  obtain  the  following  expression  for  the  general  term  in  the 
Taylor  expansion  of  a  two-centre  function: 


3/2 

n!  q,£i,£2,mi,m2,L],L2,Mi,M2,h,M 


1/2 

(».*.0OUO)  (L.L,00|L0) 


(£l,00|  A0)(£i£2mim2|ilm)  (L1L2M1M2  |  LM)  (£LmM|  Am).)j^j  ^  (ri,r2).  (4.12 
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Altt'rnatc  expressions  lor  the  expansion  about  two  eentres  are 
possible.  Wc  {)resent  the  formula  for  the  general  term  of  the  Taylor 
series  for  one  such  modification.  In  particular,  when  the  vectorial 
difference  I'l  -  r2  =  R,  a  vector  drawn  between  two  points  of  reference, 
is  a  constant  quantity,  then  the  expansion  can  lie  considered  in  terms 
of  displacements  about  the  two  ends  of  the  vector  R.  The  formula 
for  the  general  term  is 


i  [a-Vi  +  b-V2  J'VhRJ  =  I  (-i'l 

q,£i  ,^2  .111!  ,i'i2 


n  +  L 


q-  q,ii2 


21,+ 1 

[  ( 2  £  iTTynrrni  r  2  a +T) 


1/2 


(1:  1  ^  2  00|  an  iu,00|  AO)  l£,£2mim2  |£in)  (£l.mM|  A|i)V^,  (  a  )  Y  ^  ^  ( I) ) 


I  1  .  1  :n 


in  which  'S  defined  by  eqn  (2.14).  The  analysis  which  leads 

to  ec|n  (4.15)  parallels  the  proceeding  analyses. 


5.  Discussion 


Many  specific  examples  of  functions  which  can  he  expanded 
usefully  in  I'aylor  scries  easily  spring  to  mind.  Several  examines 
which  can  be  derived  from  the  Yukawa  potential  arc  considered  in  a 
companion  paper  (McKinley  and  Schmidt,  198  a). 

Wc  conclude  this  paper  by  showing  the  limiting  forms  for  the 
Taylor  scries  when  the  function  to  be  expanded  is  scalar. 


1 
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Thc  scalar  function  is  expressed  as 


CCr)  =  Voo(r3 [/^uPCr)] 

=  Hr). 


(S.l) 


The  I'ourier  transform  of  this  function  is  i'o  o  (.!<)  f  (.3;)  with  f(kj 
given  by 


f(.k)  =  (.4  II  I  dr  r^f  (r)  j  0  (kr)  , 
^  0 


tS.21 


In  terms  of  this  form  of  scalar  function,  the  general  term  for  tlie 
Tavlor  series  is 


i,  u.VJ-'c.Cr)  .  (4,0^^^(c"/n!)  I  (' i  )  V,  ,( ? ) 

il ,  m ,  L ,  M 


M/(2b+l}  (2£+l)  (LicOO  I  00)  (L£Mm|  00)  l„L(r) 


(r)i-  .  fc) 


,,  >3/2,  n,  ,>  V  ,  .>n+£.  ■* ''  ‘  £,-m'- ‘ ‘  £m'  ,  , 

=  (4ii)  (.c  /n!)  I  C-j)  - l„o  Cr) 


£  ,m 
,n+£. 


2£  +  l 


n£ 


=  /^(c"/n!  )^(-i)"  Aj^^Pj^(r*c)  (r) 


n£ 


Coulomb  potential  in  a  Legendre  polnomial  series: 

I  r+c  I  ‘\Cr*c)  . 


The  general  differential  form  also  assumes  a  simpler  ap]icarance 
in  the  scalar  limit.  We  need  only  consider  the  quantity  in 


'  21  - 


ciiii  tS.3).  I  rom  ttio  clofinition  of  wt'  liave 


elk  k''‘'“f(kjj.(,kr) 


cl 

rcTi' 


i  (cl/Jr)"’ Vl-(r) 


(5.5) 


lollowing  I'odci ,  K;iy  and  Silverstone  (1970),  wc'  can  show  that 


ef  d 

?Tr 


'r'\d/dr)"''^rP--(r) 


y  (-l)‘'(ii+q)  ! 

q  =  0  (^-mT!  (2ci‘)  !  ! 


n  -  Cl  d 

T  dr 


x(cl/Jrj'^''''^l'(r) 


I  5.0) 


TIuis  , 


'..an  = 


n  V 


1  .an^£J/2  f  (-l)‘'(£ni)! 

>r,  '  '  Jo  tF^nrniTrr 


-M 


n  -  Cl  d 

■  r  J7 


(cl/dn''''*‘^r(r). 


1 5 .  n 


lor  example,  the  second  order  term  in  the  laylor  series  assumes 
a  s  I  liip  1  r  to  nil :  \  i  , 


1  2  dl  ^J- 
•M  r  a?  j;. 


=  1’  -  (  r -c  ) 
,•)  ■ 


1  ell- 
r  cl  r 


dji 
d  1-  n 


(  5  .  .5  ) 


This  form  is  particularly  useful  foi  .i  number  ol  applications  whicli 
involve  scalar  functions. 
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